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Abstract:  One  of  the  important  applications  of  underwater 
communication  is  speech  transmission  between  two  divers  or 
between  divers  and  ship  or  submarine.  This  paper  describes  a 
project  designed  to  investigate  and  demonstrate  underwater 
communication  system  in  Persian  Gulf  for  speech  transmission  in 
a  real  channel.  At  first,  transmitter  is  designed,  then  channel  with 
real  data  is  simulated  by  neural  network  and  at  last  receiver  is 
designed.  Transmitted  data  is  speech  signal  that  for  more  safety 
in  transmission  and  low  frequency  bandwidth,  a  cryptography 
algorithm  and  speech  coding  algorithm  is  applied  in  transmitter. 
Quadrature  phase  shift  keying  (QPSK)  signaling  is  employed  to 
make  efficient  use  of  the  available  channel  bandwidth.  In  the 
receiver,  linear  equalizer  and  decision  feedback  equalizer  (DFE) 
are  tested  and  the  best  of  them  is  applied.  Also  ray  tracing 
method  is  used  for  simulation  of  sound  waves  propagation  in 
Persian  Gulf  underwater  communication  channel. 

L  Introduction 

The  underwater  acoustic  communication  channels,  however, 
have  limited  bandwidth,  and  often  cause  signal  dispersion  in 
time  and  frequency  [1].  In  this  paper  we  have  designed 
transmitter  and  receiver  that  transmit  speech  signal  in  the  safe 
mode  and  low  bandwidth  with  low  bit  error  rate  (BER)  in 
Persian  Gulf  underwater  communication  channel. 

In  the  first  section,  transmitter  is  designed.  Because  of 
strong  frequency  attenuation,  channel  bandwidth  is  limited. 
Therefore  in  transmitter,  we  have  used  LPC-10  (Linear 
Prediction  Coder- 10)  algorithm  to  compressing  speech  signal. 
After  compressing,  for  safety  in  transmission,  RC5  (Rivest 
Cipher)  cryptography  algorithm  has  used  to  encrypt  data.  Then, 
for  reduction  of  BER,  channel  coding  is  used.  In  the  last 
section  of  transmission,  data  is  modulated. 

In  third  section,  Persian  Gulf  underwater  channel  is 
simulated.  In  channel  modeling,  speed  and  depth  profiles  of 
Persian  Gulf  have  used  that  have  measured  in  358  zone  in 
Persian  Gulf  by  Canadian  in  2001. 

In  fourth  section,  linear  equalizer  and  DFE  equalizer  are 
designed. 

In  fifth  section,  transmitter,  Persian  Gulf  communication 
channel  and  receiver  are  simulated  for  real  condition  and  in 
the  last  section,  the  results  are  presented. 

II.  Transmitter 

In  this  section,  transmitter  has  designed.  Fig.  1,  shows 
transmitter  block  diagram  [2].  In  the  continuation,  we  will 
describe  any  block. 


Figure  1 .  Transmitter  bloek  diagram 


A.  Speech  signal  processing 

In  the  first  block  of  transmitter,  speech  signal  is  compressed 
with  LPC-10  algorithm.  Sampling  frequency  is  8KHz  and  for 
any  sample  8  bit  is  appropriated.  Therefore  sampling  bit  rate  is 
64  kbps.  LPC-10  algorithm  reduce  bit  rate  to  2.4  kbps, 
consequently  small  bandwidth  is  needed  and  frequency 
attenuation  is  reduced  [3,  4]. 

B.  Cryptography 

In  the  second  block,  data  is  encrypted  with  RC5  algorithm. 
RC5  is  a  symmetric  block  cipher  designed  to  be  suitable  for 
both  software  and  hardware  implementation.  It  is  a 
parameterized  algorithm,  with  a  variable  block  size,  a  variable 
number  of  rounds  and  a  variable-length  key. 

This  provides  the  opportunity  for  great  flexibility  in  both 
performance  characteristics  and  the  level  of  security.  In  this 
paper,  block  size  has  selected  32,  number  of  rounds  are  16  and 
key  length  is  10  [5]. 

C.  Channel  coding 

This  block  is  needed  for  error  correction  in  channel.  BCH 
(Bose,  ray-Chaudhuri,  Hocquenghem)  coding  has  been 
selected  for  channel  coding.  The  inputs  of  applied  BCH  code 
are  11  bits  that  convert  to  31  bits  and  are  transmitted  in 
channel.  This  code  can  correct  5  or  fewer  random  errors  in 
receiver  [6]. 

D.  Modulator 

Modulator  utilize  quadrature  phase  shift  keyed  (QPSK) 
modulation.  The  encoder  accepts  the  sequence  of  the  input 
binary  data.  It  has  two  outputs;  in-phase,  I,  and  quadrature,  Q. 
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For  each  distinct  pair  of  input  binary  data  a  unique 
combination  of  /  =  j+ljand^  —  assigned.  We  can 

use  a  complex  notation  D  =  I  +  jQ  to  describe  these  four  points 
separated  in  the  complex  plane  by  45% 135% 22 5°  or  315°  . 
The  encoded  data  stream  of  I  and  Q  is  then  used  to  modulate  a 
sequence  of  S  pulses  Transmitted  every  signaling  period  T. 
To  limit  their  bandwidth  such  modulated  sequences  are  then 
low-pass  filtered  by  LPFs.  The  same  low  pass  filters  are 
applied  at  the  receiver.  The  in-phase  and  quadrature  signals  at 
the  output  of  the  transmitting  low  pass  filter  are: 


Time(msec) 


Figure  2.  Impulse  response  of  the  raised  eosine  filter 


'  (1) 

Ssit)=XQ^hAt-iT) 

i 

where  {t)  is  the  impulse  response  of  the  low  pass  filter. 
The  filtered  signals  are  then  multiplied  by  a  carrier  frequency, 
added  and  transmitted  through  the  Persian  Gulf  underwater 
acoustic  channel. 

In  order  to  reduce  inter-symbol  interference  (ISI)  we  use 
modified  raised-cosine  filters  (LPFs)  with  a  roll-off 
factor  P  =  \  and  impulse  response  {t)  .  The  transfer 
function  of  the  baseband  channel  has  the  form: 

(2) 

0  ;  elsewhere 

where  T  =  1/1200  s.  The  transfer  function  of  both  the 
transmitter  filter  and  receiver  filter  is  then  f ),  and  the 
impulse  response  of  the  assumed  nonfading  part  of  the 
baseband  channel  is  thus: 


To  reduce  sidelobe  levels,  the  impulse  response  of  the  filter 
is  modified  by  multiplying  it  with  a  Hamming  window  given 
as: 

(0=  0.54  +  0.46  cos(  (5) 

Fig.  3  shows  Hamming  window  in  one  period.  The  period 
of  Hamming  window  is  2T  that  T  is  the  period  of  each  symbol. 

This  window  has  99.96%  of  its  energy  in  the  main  lobe, 
with  sidelobes  of  over  20  dB  down  from  the  peak,  as  shown  in 
Fig.  4  [7,  8]. 


Figure  3.  Hamming  window  in  one  period 


x{t)^ 


(3) 


Thus  the  impulse  response  of  each  filter  is  obtained  through 
the  inverse  Fourier  transform  of  Eq.  (2),  which  is: 
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Fig.  2  shows  the  impulse  response  of  the  raised  cosine  filter. 


Frequency  (Hz) 

Figure  4.  Impulse  response  of  LPF  with  and  without 
Hamming  window 
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Thus  Impulse  response  of  filters  is  corrected: 

(0  ~  (0*^//  (0  ~  (0*^//  (0  (^) 

We  use  this  modified  raised  cosine  filter  in  the  transmitter 
and  receiver.  LPF  outputs  in  transmitter  are: 


Structure  of  this  neural  network  has  shown  in  Fig.  7.  This 
network  has  two  hidden  layer  with  10  neurons  in  each  layer. 

Too,  other  neural  network  has  been  designed  for  sound 
celerity  estimation.  The  inputs  of  this  neural  network  are 
longitude,  latitude  and  depth  and  its  output  is  sound  celerity. 
With  this  neural  network  we  can  estimate  sound  celerity  in 
every  point.  Neural  network  structure  is  shown  in  Fig.  8. 


(7) 

SQ{t)=Y.QihHit-iT) 

i 


For  example  if  the  input  of  LPF  may  be  sequence 
101 110001 1  then  output  has  been  shown  in  Fig.  5.  The  output 
signal  of  modulator  has  been  shown  in  Fig.  6. 

III.  Channel  simulation 

There  are  many  methods  for  acoustic  wave  propagation  in 
underwater  channel  [9].  In  this  paper  we  have  used  ray  method 
for  acoustic  wave  propagation  in  Persian  Gulf  underwater 
channel  [10,  11]. 

The  salinity,  depth  and  sound  speed  have  been  measured  by 
Canadian  in  358  points  with  different  longitudes  and  latitudes 
in  2001.  This  information  for  other  points  has  been  estimated 
by  neural  network.  At  the  first,  we  have  designed  a  neural 
network  that  its  inputs  are  longitudes  and  latitudes  and  its 
output  is  depth.  Therefore  with  this  neural  network  we  can 
estimate  depth  in  every  point. 


Figure  5.  Output  of  LPF  for  sequence  1011 10001 1 


IV.  Receiver  design 

The  digital  acoustic  communication  is  limited  by  severe 
inter-symbol  interference  (ISI)  associated  with  long  range, 
shallow  underwater  channels.  The  ISI  is  caused  by  multi-path 
propagation  resulting  from  surface  and  bottom  reflections.  To 
overcome  the  effects  of  intersymbol  interference,  an  adaptive 
equalization  method  employing  a  Mean  Square  Error  (MSE) 
criterion  is  introduced.  Computer  simulations  are  carried  out 
to  verify  the  effectiveness  of  the  equalization  scheme  for  high 
data  rate  communication. 

In  the  selection  of  an  equalizer  we  consider  convergence 
rate,  processing  complexity  and  the  ability  to  track  changes  in 
the  channel  characteristics. 


Input  Hidden  Hidden  Output 
Layer  Layerl  Layer2  Layer 

Figure  7.  Designed  neural  network  for  depth  estimation 


0  0  0  0 


Input  Hidden  Hidden  Output 

Layer  Layerl  Layer2  Layer 


Figure  6.  Output  signal  of  modulator 


Figure  8.  Designed  neural  network  for  sound  eelerity  estimation 
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There  are  a  variety  of  algorithms  which  have  been 
investigated  for  use  in  underwater  acoustic  communication. 

In  this  paper  we  have  tested  linear  and  DFE  equalizers  for 
speech  signal  transmitting  in  Persian  Gulf  underwater  channel. 


New  set  of  coefficients  is  obtained  iteratively  as: 

W|  [{k  +  l)r]  =  W|  {kT)  +  ^e{kT)y*  {kT  -  ir)  (10) 


A.  Linear  equalizer 

A  symbol-spaced  linear  equalizer  consists  of  a  tapped  delay 
line  that  stores  samples  from  the  input  signal.  Once  per 
symbol  period,  the  equalizer  outputs  a  weighted  sum  of  the 
values  in  the  delay  line  and  updates  the  weights  to  prepare  for 
the  next  symbol  period.  This  class  of  equalizer  is  called 
"symbol-spaced"  because  the  sample  rates  of  the  input  and 
output  are  equal. 

Fig.  9  is  a  schematic  of  a  symbol-spaced  linear  equalizer 
with  L  weights,  where  the  symbol  period  is  T. 

Impulse  response  of  linear  equalizer  is: 

+  (8) 

We  should  determine  {H;}^_^to  minimize  error  probability. 

The  criterion  used  in  the  optimization  of  the  equalizer 
coefficients  is  the  minimization  of  the  mean  square  error 
(MSB),  that  is  [12]: 

min  {E\e{kT  f  |)= 

(T  D 

min  ^E\z^{kT)-z{kTf^ 

Where  Zd(kT)  and  z(kT)  are  the  actual  transmitted  data  and 
the  equalizer  output  at  t=  kT,  respectively.  For  an  ideal 
channel,  d(kT)  =  z(kT).  During  the  equalizer  training  period,  a 
block  of  data  known  to  the  receiver  can  be  transmitted  and 
used  to  adjust  the  equalizer  coefficients.  After  the  training 
period  a  receiver  does  not  know  the  actual  transmitted  data  so 
that  the  adjustment  of  the  equalizer  coefficients  (tap  gains)  is 
performed  using  the  decisions  at  the  output  of  the  detector 
(decision-directed).  This  is  symbolized  by  switch  S  in  Fig.  9. 
The  minimization  of  the  MSB  is  accomplished  recursively  by 
the  use  of  the  stochastic  gradient  algorithm,  called  the  least 
mean  square  (BMS)  algorithm. 


where  A  is  a  positive  adaptation  constant  (step  size), 
e(kn)  =Zd(kT)  -  z(kT)  is  the  error  signal  and  *  denotes  the 
complex  conjugate. 

The  step-size  parameter  A  controls  the  rate  of  adaptation  of 
the  equalizer  and  stability  of  the  stochastic  gradient  algorithm. 
Barger  A  leads  to  more  rapid  convergence,  but  also  leads  to 
larger  fluctuations  in  the  equalizer  coefficients  during  steady 
state  operation.  The  selection  of  A  is  therefore  a  compromise 
between  rapid  convergence  and  small  variance  of  self  noise. 

B.  Decision  feedback  equalizer 

A  decision-feedback  equalizer  is  a  nonlinear  equalizer  that 
contains  a  forward  filter  and  a  feedback  filter.  The  forward 
filter  is  fractionally  spaced  with  spacing  T/2;  while  the 
feedback  filter  contains  a  tapped  delay  line  whose  inputs  are 
the  decisions  made  on  the  equalized  signal.  The  purpose  of  a 
DFE  is  to  cancel  intersymbol  interference  while  minimizing 
noise  enhancement.  By  contrast,  noise  enhancement  is  a 
typical  problem  with  the  linear  equalizers  described  earlier. 
Fig.  10  shows  structure  of  decision  feedback  equalizer  (DFE). 

Similar  to  linear  equalizer,  the  criterion  used  in  the 
optimization  of  the  equalizer  coefficients  is  the  minimization 
of  the  mean  square  error  (MSB)  by  use  of  the  least  mean 
square  (BMS)  algorithm.  The  output  of  DFE  at  t=kT  is  given 
as  [13,  14,  15]: 


N-\ 

z{kT)  =  Y,  c^{kT  )y{kT  -  nT  !  l)- 


n  =  0 
M 

^b^{kT)dikT  -m) 

m=\ 


(11) 


And  error  is: 

e{kT)  =  d{kT)  —  z{kT)  (12) 


Figure  10.  Structure  of  decision  feedback  equalizer 
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New  set  of  coefficients  is  obtained  iteratively  as: 


c.[(k+i)r]= 

c„  {kl)  +  ^e{kT)y  {kT-nT/2) 

b.[{k  +  l)T]  = 

{kT)  +  ^e{kT)d*{kT  -niT) 


(13) 


(14) 


VI.  Simulation  results 

In  this  section,  we  simulate  Persian  Gulf  underwater 
channel  in  September  2001  and  both  linear  and  DFE 
equalizers  are  tested  in  receiver.  Table  I  shows  the  parameters 
of  the  channel  and  system. 


TABLE  I 

SIMULATION  PARAMETERS 


Range 

10  km 

Water  depth 

120  m 

Transmitter  depth 

15  m 

Receiver  depth 

30  m 

Carrier  frequency 

10  kHz 

Bandwidth 

2  kHz 

Fig.  11  shows  the  sound  speed  profile  in  Persian  Gulf  in 
September  2001  and  Fig.  12  shows  directions  of  sound  wave 
propagation  in  Persian  Gulf  underwater  channel.  Fig.  13 
shows  the  length  of  any  wave  propagation  path  vs.  meter  and 
Fig.  14  shows  the  attenuation  of  any  path.  The  longer  paths 
have  more  attenuation.  The  path  that  is  longer,  its  attenuation 
is  more.  By  dividing  rang  to  sound  speed  in  water,  we  can 
calculate  the  delay  of  each  path. 


Depth  profile 


■jupBTr:  -  lim 
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Range(m) 

Figure  12.  Direction  of  sound  wave  propagation  in  Persian  Gulf  underwater 

channel 


1  2  3  4  5  6  7  8  9  10  11  12  13 

Path  number 

Figure  13.  Range  of  paths 


Speed(m/s) 


Figure  1 1 .  Speed  profile  in  September  2001  in  Persian  Gulf  at  longitude  59^ 
and  latitude  24® 


Path  number 


Figure  14.  Attenuation  of  paths 

Structure  of  the  receiver  has  been  shown  in  Fig.  15.  We  can 
replace  adaptive  equalizer  block  with  linear  or  nonlinear  (DFE) 
equalizer.  The  output  of  channel  r(t)  and  outputs  of  EPF  in 
receiver  are  shown  in  Fig.  16,  Fig.  17  and  Fig.  18. 
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Figure  15.  Structure  of  the  receiver 


1.15  1.2  1.25  1.3  1.35 


Time(sec) 

Figure  16.  Output  of  channel 


Time(sec) 


Figure  17.  Output  of  LPF  in  receiver 


Time(sec) 

Figure  18.  Output  of  LPF  in  receiver 

With  respect  to  maximum  delay  in  channel,  linear  equalizer 
has  46  weights  and  A=  10'"^.  DFE  has  92  weights  in  forward 
filter  and  46  weights  in  backward  filter. 

2500  symbol  is  applied  for  equalizers  training.  As  can  be 
seen  in  Fig.  19,  for  convergence,  the  equalizers  requires  over 
100  iterations. 


The  output  of  linear  and  DFE  equalizer  in  SNR=0  dB  have 
been  shown  in  Fig.  20  and  Fig.  21. 

In  Fig.  22  and  Fig.  23,  we  have  shown  speech  signal  in 
transmitter  and  receiver. 


Figure  19.  BER  vs.  number  of  iterations 


-5  0  5 

Inphase 


Figure  20.  Output  of  linear  equalizer  in  SNR=0dB 
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Inphase 

Figure  21.  Output  of  DFE  in  SNR=0dB 


Figure  22.  Speech  signal  in  the  transmitter 
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Figure  23.  Speech  signal  in  the  receiver 


VII.  Results 

In  this  paper,  we  designed  an  underwater  communication 
system  in  Persian  Gulf.  In  transmitter,  data  are  compressed, 
encrypted  and  modulated  for  sending  in  underwater  channel. 
Then  we  simulated  wave  propagation  in  Persian  Gulf 
underwater  channel  by  using  neural  network  and  ray  method. 
In  receiver,  we  tested  linear  and  decision  feedback  equalizer. 
Results  show  that  DFE  is  suitable  for  Persian  Gulf  underwater 
channel.  Also  data  are  decrypted  and  decoded. 
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